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Simplex Spline Bases on the Powell-Sabin 12-Split: Part I 
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(joint work with Elaine Cohen, Richard Riesenfeld) 

Abstract. We review the construction and a few properties of the S- 
basis, a simplex spline basis for the quadratic splines on the Powell- 
Sabin 12-split. 

Piecewise polynomials or splines defined over triangulations form an indispensable 
tool in the sciences, with applications ranging from scattered data fitting to finding 
numerical solutions to partial differential equations. In applications like geometric 
modelling and solving PDEs by isogeometric methods one often desires a low 
degree spline with C^, or smoothness. For a general triangulation, it is 
known that the minimal degree of a triangular C element is 4r -|- 1, e.g., degrees 
5, 9, 13 for the classes G^, G^ or G^. To obtain smooth splines of lower degree 
one can split each triangle in the triangulation into several subtriangles. One such 
split that we consider here is the Powell-Sabin 12-split of a triangle. 



The 12-split with numbering of vertices. 

Once a space is chosen one determines its dimension. The spaces and 

iSKA) of G^ quadratics and G^ quintics on the 12-split A of a single triangle have 
dimension 12 and 39, respectively. Over a general triangulation T of a polygonal 
domain we can 12-split each triangle in T to obtain a triangulation Tn- The 
dimensions of the corresponding G^ quadratic and G^ quintic spaces (the latter 
with G^ supersmoothness at the vertices and the interior edges of each macro 
triangle) are 3|V| -I- \£\ and 10|V| -I- 3|£|, respectively, where |V| and \£\ are the 
number of vertices and edges in T. Moreover, in addition to giving G^ and G^ 
spaces on any triangulation these spaces are suitable for multiresolution analysis, 
see for example [5]. 

To compute with these spaces one needs a suitable basis. In the univariate 
case the B-spline basis is an obvious choice. In this talk we consider a bivariate 
generalization known as simplex splines. We review the construction and a few 
properties shown in [1] of the S-basis consisting of G^ quadratic simplex splines in 
52 (A). We also introduce some concepts needed in Part II of this talk given by 
Georg Muntingh. 

A short background on simplex splines. Let K = {ui • • • u^-i-s+i} C K® 
be a finite multiset. Consider a simplex a = [ui, ..., C together 
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with a projection tt : a —> K® satisfying TT{vi) = Vi. We define the simplex spline 
B[K]{x) = vold(cr n 7r“i(a;))/vold+s (cr). For instance, three knots in define a 
linear B-spline, four knots in define a quadratic B-spline, and four knots in 
define a linear bivariate simplex spline: 



Simplex splines have all the usual properties of univariate B-splines. This in¬ 
cludes continuity which can be controlled locally, a recurrence relation, and dif¬ 
ferentiation and knot insertion formulas. The support of a simplex spline is the 
convex hull of its knots, and in the collection of knotlines is obtained by con¬ 
necting each knot to all other knots (the complete graph). A simplex spline with 
d + 3 knots in has d — m + 1 continuous derivatives across a knot line containing 
m knots counting multiplicites. 

Simplex splines on the 12-split. Since the knotlines form a complete graph 
the simplex splines are natural candidates for a C"' basis on this split. A simplex 
spline on the 12-split will have a knotset of the form K = where 

Di,... , 1)10 are the vertices numbered as above, and rui > 0 is the multiplicity of 
Vi, i.e., the number of repetitions of Vi in the multiset. A convenient scaling is the 
(area normalized) simplex spline Q[K] : —> K, recursively defined by 


Q[K]{x) 


0 

1 / \ area(Z\) 

[JS^)(®)area([JC]) 


if area([iir]) = 0, 
if area([iir]) ^ 0 and |lir| = 3, 
if area([iir]) ^ 0 and |iir| > 3, 


with X = PiVi -!-■••+ PiqViq, j3l^ --I- /3io = 1, and /3i = 0 whenever mi = 0. 

By Theorem 4 in [3] this definition is independent of the choice of the jdj. 
Whenever m-j = ms = mg = mig = 0, we use the graphical notation 



:= Q\vlvivlv[v^v^] 


Q[ 


B-splines on the boundary. It is useful for the simplex splines to restrict 
to consecutive univariate B-splines on the boundary. For example, on [vi,V 2 ] the 

quadratic simplex splines I restrict to: 
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form a group S 3 that acts on the simplex splines by permuting knots. We write 
[B]s, := {Q[a{K)] : Q[K] & B, a e S3} 
for the set of simplex splines related to S by a symmetry in S3. Let 


C4 := 

C7 := 


Cl + C2 
2 ^ 
C4 + Cg 


C 5 := 
C8 := 


C 2 + C 3 
2 ’ 
C4 + C5 


Cg := 
Cg := 


Cl + C3 
2 ’ 
Cg + Cg 


cio := 


Cl + C2 + C3 


Via the identification Ci O Vi with the vertices of A, the group S3 acts on poly¬ 
nomials in Cl,, cio and simplex splines, or combinations of these, e.g.. 



The quadratic S-basis. It is given by 



and is the unique simplex spline basis for 52 (A) with local linear independence. 
Moreover, it is symmetric, reduces to B-splines on the boundary, can be computed 
by a pyramidal scheme, and has Bezier-like smoothness conditions across adjacent 
macro triangles. Furthermore, it has a barycentric Marsden identity 



which yields polynomial reproduction, explicit dual functionals and a simple quasi- 
interpolant. These show that the S-basis is stable independently of the geometry, 
which implies an bound on the distance between a spline and its control surface. 
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